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DIGEST 


Condensation  processes  are  important  in  the  thermal  dissemination 
of  CIV  agents.  Heterogeneous  3nd  homogeneous  nucleatior  was  stud’rd  tc 
determine  the  important  factors  controlling  condensing  systems. 

F. ffects  of  salt  nuclei  on  the  parLicle  size  distribution  of  the 
disseminated  aerosol  were  studied  and  special  pyrotechnic  systems  which 
were  salt  nuclei-free  were  investigated. 

Aerosols  composed  of  two-  to  five- micron-diameter  particles  with  a 
high  degree  of  particle  size  homogereity  make  the  optimum  aerosol  for 
lung  retent ion  and  maximum  transparency.  The  production  of  such  an 
aeiosol  was  shown  to  be  feasible  bv  the  technique  of  incorporating  non¬ 
volatile  "giant  nuclei"  material  in  the  pyrotechnic  mix.  These  giant 
nuclei  when  disseminated  with  the  agent  vapor  aited  as  preferential 
condensation  sites  and  as  small  particle  scavengers  by  coagulation. 

The  secondary  process  of  coagulation  was  shown  to  be  important  in  re¬ 
moving  the  highly  visible  submicron  particles.  A  pyrotechnic  dissem¬ 
ination  system  was  suggested  to  produce  uniform,  low  visibility  aerosols 
which  included  coagulation  of  the  small  particles  to  be  disseminated  on 
giant  nuclei  at  elevated  temperature  and  high  concentrations. 

Homogeneous  nucleation  always  occurs  in  condensing  systems  of  hi  gh 
vapor  concentrations  even  in  the  presence  of  foreign  nuclei.  Critical 
supersaturation  ratios  of  a  number  of  compounds  were  measured  by  a 
newly  developed  experimental  method.  It  was  demonstrated  that  the 
classical  Be  eke r - Doe r i ng  theory  is  inadequate  for  the  preparation  of 
nucleation  models.  A  modified  model  is  presented  showing  the  inter¬ 
relationships  among  physical  properties,  chemical  structure,  and 
temperature . 
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I  INTRODUCTION 


Pyrotechnic  dissemination  ol  liquid  or  solid  agent  is  accomplished 
by  mixing  (lie  agent  with  the  pyrotechnic  materia).  The  phase  changes  in 
the  process  include  volatilization  to)  lowed  by  cnndt'iissi  of  th._  agent - 
The  manner  in  which  condensation  occurs  establishes  t  lit  particle  size 
distribution  of  the  disseminated  material.  11  the  condensation  processes 
in  pyrotechnic  dissemination  could  bt  isolated  and  studied  separately, 
then  the  agent  paiticle  size  and  homogeneity  could  bo  controlled. 

Condensation  involves  nucleation,  growth,  and  coagulation  (agglom¬ 
eration).  Se I f-nuc 1  eel i on  (homogeneous  nucleation)  always  occurs  at  some 
critical  supersaturation  of  a  condensable  vapor,  but  foreign  nuclei,  if 
they  are  active,  will  initiate  condensation  at  lower  supers.ilurat  ion  values. 
The  conditions  of  agent  aerosol  formation  in  pyrotechnic  dissemination  are 
highly  specific,  and  no  theoretical  analysis  of  vapor  condensation  under 
these  conditions  has  ever  been  attempted.  lhe  size  distribution  of  the 
agent  aerosol  involves  the  si  mu  I t  ..neous  kinetics  of  nucleation,  growth, 
and  agglomeration. 

Condensation  in  rapidly  quenched  jits  i »  dependent  on  the  number  ana 
kind  of  nuclei  present.  There  is  considerable  experimental  evidence  that 
under  rapid  quench  conditions,  in  the  presence  ol  foreign  nuclei,  homo¬ 
geneous  and  heterogeneous  nucleation  occur  simultaneously.  Therefore,  in 
order  to  understand  condensation  in  quenched  jets  it  is  necessary  to  study 
botli  heterogeneous  and  homogeneous  nucleation  processes.  He  have  made  ex¬ 
perimental  studies  ol  both  types  of  nucleation.  These  have  included 
(1)  an  investigation  ol  the  aerosol  properties  of  material  disseminated 
from  a  pyrotechnic  device  and  (2)  measurement  of  the  homogeneous  nuclea¬ 
tion  oi  a  number  of  simple  compounds  to  understand  more  fully  the  importance 
ol  homogeneous  nucleation  in  actual  systems. 
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II  NUCLEATION  IN  PYROTECHNIC  DISSEMINATION 


A.  Isomorphic  Nuclei  Kf feels 

The  effect  o  1  foreign  mu  In  <>:i  the  particle  size  distributions  of 
solid  agent  si.nulunls  vaporized  from  pyrotechnic  devices  was  studied  using 
u  iiunihci  of  different  blends  of  pyrotechnics  and  agent  simulants.  These 
mixtures  were  burned  in  unpacked  piles  and  the  disseminated  simulant  '.as 
analyzed  for  size  distribution  using  a  six-stage  impactor,  the  Andersen 
Sdipler.*  Benzoic  acid  (BA)  phthalic  anhydride  (PA)  and  1 -methyl  - 
an noan t h raqu i none  (MAA)  wer>  the  principal  simulants  used.  Potassium 
chlorate  (KClOj)  and  ammonium  piciato  (APic)  were  the  oxidizers  used. 

Inorganic  s  ,  1 1 s  ( 5%  by  weight!  were  added  to  the  AP i c • s i mu  1  an t  mixes 
to  provide  suit  nuclei  upon  burning.  The  effect  o f  salt  nuclei  can  be 
noted  in  Table  I.  As  in  the  case  of  nucleation  in  condensing  jets, 
large  numbers  of  salt  nuclei  tend  to  dec  tease  the  particle  size  and  the 
particle  size  range  of  the  recondensed  material.  Also,  as  observed  in 
the  jet  nucleation,  the  salts  which  arp  isomorphic  appear  to  produce  the 
greatest  effets.  The  foreign  nuclei  concentration  in  these  systems  is 
probably  critical  and  may  be  a  quantitative  function  of  nucleation 

efficiency. 

The  effectiveness  of  selected  isomorphic  nuclei  upon  condensation 
is  demonstrated  by  the  comparative  particle  sizes,  but  the  final  aerosol 
propertifs  depend  on  the  relative  rates  of  nucleation,  growth  from  the 
vapor,  and  agglomeration  of  particles  already  formed.  Thus,  it  would 
appear  to  be  possible  to  control  size  of  the  particles  in  a  thermal 
d i s semi n a t  ion  jet. 


i 


manufactured  by  the  Andcraen  Sampler  O.  , 
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*r>  Nor  re.  oitf  dlllnlijUo.l. 

B.  Thermal  Dissemination  with  Giaru  Nuclei 

The  possibility  of  controlling  the  particle  size  of  O  agent  dis¬ 
seminated  by  a  pyiotechnic  device  was  investigated  experimentally.  In 
order  to  control  particle  size,  giant  nuclei  were  added  and  mixed  into 
the  py  re.:  ccl.r.  ic  -  agent  blend  The  giant  nor  lei  rr.alci  i  «  1  had  to  have  very 
definite  physical  properties: 

(1)  particle  sizes  had  to  be  in  the  !-  to  3-rr.icron- 
diameter  range, 

(2)  the  material  cooid  n<t  volatilize  appreciably  at 
pyrotechnic  burning  temperatures;  and 

(3)  the  material  had  to  be  a  preferential  nucleator 
(small  li.teiiuciai  tension  with  agent). 

The  materials  investigated  were  silica,  magnesium  oxide,  titanium  dioxide 
ferric  oxide,  and  alumina. 

The  pyrotechnic  mix  could  not  be  a  source  of  effective  nuclei.  Tins 
last  condition  was  the  most  difficult  one  to  attain.  Ammonium  piciate,  a 
ir.onopropc  1  1  an t  which  contains  no  melui  atom  and  produces  no  salt  nuclei, 
was  the  lirst  fuel- oxidizer  used. 
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The  dependence  of  particle  size  of  condensed  organic  aerosol  particles 
on  the  presence  of  giant  nuclei,  which  were  added  to  the  bulk  mat. 'rial 
before  vaporization,  was  demonstrated  in  a  simple  n onp y rotec hn i c  volati¬ 
lization.  In  an  attempt  to  make  a  comparative  study  of  the  particle  size 
distribution  of  organic  simulants  their.ially  disseminated,  with  and  without 
selected  giant  nuclei  added  to  the  pyrotechnic  mix,  a  number  cf  mixes  were 
blended  and  pressed.  All  mixes  contained  APic  as  the  monopropellant. 

The  mixes  were  made  with  APir  (50  and  7 O'T-  by  weight)  using  benzoic  acid, 
phthalic  anhydride,  and  1  -  me t hy 1  ami noan t h requ i none  as  agent  simulants. 
Additional  mixes  were  made  with  the  same  ingredients  plus  giant  nuclei 
material  to  the  extent  of  1 0‘7  by  weight  of  'he  simulant.  The  giant 
nuclei  materials  were  1-  to  4 . 5 - m i c r on - d l ame t e r  silica  and  a  pigment 
grade  titanium  dioxide.  The  mixes  (i.e.,  nuclei,  propellant,  and  agent- 
si  mi  1  e  :u  )  were  blended  in  a  shaker  for  several  hours  and  then  pressed. 

A  50:50  mixture  of  APic  and  BA  and  the  same  mixture  plus  55c  1-  to  4-micron- 
dia meter  silica  were  fired  in  s  e  p  rate  experiments  in  a  one- cub ■ c- foot 
meter  box.  The  smoke  from  each  mixture  was  collected  in  the  Andeisen 
Sampler  plus  Millipore  filter  train  The  relative  amounts  of  each  mix¬ 
ture  collected  ori  the  backup  filters  showed  that  there  was  several  times 
as  much  finely  divided  material  from  the  50  SO  mixture  as  from  the 
mixture  containing  silica.  However,  in  both  cases  there  was  black  smoke, 
presumably  from  the  inefficient  combustion  of  the  APic. 

Pyrotechnics  based  on  potassium  chlorate  will  supply  large  numbers 
of  finely  divided  (0.1  to  0. 01  -  mi c r on - di ame ter )  KCi  nuclei.  A  second 
monopropellant  which  is  salt,  nuclei-  free  is  t  r  l  ami  noguan  i  d  l  n  e  cyanoformate 
(TAGCY)  This  compound  which  was  developed  at  SHI  is  far  superior  to 
APic  as  an  effective  monopropellant.  TAGCY  mixtures  containing  1  and  2" 
APic  plus  3051  agent  simulant  were  successful  pyrotechnics.  Ten  gram 
grains  of  these  mixes  were  pressed  at  2000  lb  into  aluminum  film  cans  of 
3  cm  in.  diameter.  Addition  of  1-  to  4  -  m  i  c  ron  -  d  i  a  .re  t  e  r  silica  to  these 
pyrotechnic  mixes  produ<ed  a  marked  increase  in  th<  mass  median  diameter 
of  the  disseminated  material  and  a  rr.ore  homogeneous  size  range.  Submicron 
particles  were  produced  during  dissemination,  however,  and  a  more  effi¬ 
cient  scavenging  of  the  smoke  was  necessary. 

hour  systems  and  four  coagulation  cimmbeis  were  evaluated  lor  in¬ 
i' tensing  the  amount  of  coagulation  in  aerosols  produced  by  pyrotechnic 
means.  Changes  in  relative  cloud  densities  and  particle  size  distribu¬ 
tions  were  the  measured  parameters.  The  four  chamber  designs  studied  were : 


(1)  empty  cylinder  for  simple  confined  coagulations,  (2)  cylinder 
with  two  8-mesh  screens  (see  Fig.  1),  (3)  cylinder  with  two  orifice 
plates,  each  containing  fourteen  1/8-in.  holes  (see  Fig.  2),  and  (4)  a 
two-stage  coagulation  device,  in  which  initial  coagulation  occurs  at  hi  gh 
temperature  and  pressure  in  u  confined  volume  and  is  followed  by  simple 
confined  coagulation  l  r.  a  second  stage. 

.411  of  these  coagulation  devices  were  scaled  for  use  with  10  g  of 
pressed  mix  contained  in  a  30  -  mm- d i ame te l  canister. 

Three  types  of  pyrotechnic  mixes  were  used.  The  first  contained 
68%  TAGCY,  2%  APic,  27*5  BA,  and  3°*  silicic  acid  of  about  2-  to  10-micron- 
diameter  particle  size  for  nucleate  seeding.  The  second  and  third  mixes 
were  based  on  KC103 -  1 ac t os e  propellant  and  had  different  coolant  and 
seeding  materials.  The  second  mix  was  25%  KClOj,  11.5%  kaolin,  19.5% 
lactose,  40%  BA,  and  4%  silica  (1-  to  4 -mi  cron  diameter  seeding  material). 
The  third  mix  was  formulated  to  check  the  possibility  of  using  a  single 
additive,  silicic  acid,  for  both  the  cooling  and  seeding  functions.  The 
formulation  was  25%  KCllh,  19.5%  lactose,  44%  BA,  ~nd  11.5%  silicic  acid. 

Experiments  with  the  TAGCY  based  mix  revealed  that  screens  are  more 
efficient  in  promoting  coagulation  than  orifice  plates  and  that,  e’en 
though  the  velocity  through  the  on  f  u  e  holes  was  about  five  limes  greater 
than  the  velocity  through  the  screen,  the  dimensions  and  larger  number  of 
turbulent  cell*  caused  more  coagulation  with  the  screens  than  with  the 
orifice  plates. 

The  two-stage  coagulator  showed  a  greater  calculated  coagulation 
rate  than  the  single-stage  system.  The  two- stage  approach  shuwed  the 
advantage  over  the  single-stage  approach  of  suppressing  homogeneous 
nucleations  in  favor  of  heterogeneous  nucleation. 

Collision  rate  calculations  foiiow,  based  on  the  kinetic  theory  of 
gases  and  on  Smo 1 uchowsk  l  ' s  aeiosol  coagulation  theory.1  The  ratio  ol 
the  calculated  collision  rates,  confined  coagulation  (two-stage)  to 
unconfined  coagulation  (single-stage),  is  about  20:1  in  both  calcula¬ 
tions.  The  absolute  values  of  the  collision  rates  based  on  the  two 
theories  give  rales  of  the  same  magnitude. 
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74  HOLE 


DEVICE  No.  I 


DEVICE  No.2 


FIG.  1  COAGULATION 
CHAMBER  WITH 
SCREENS 


DEVICE  No.  2 


FIG.  2  COAGULATION  CHAMBER 
WITH  ORIFICE  PLATES 


The  calculations  are  based  on  the  following  assumptions: 

1.  For  every  10-g  charge  there  is  0.4  g  of  silica  seeding 
material.  Using  an  average  size  of  2  microns  in 
diameter,  the  total  number  of  seeding  particles  is 

5  x  1 0 1  In  the  confined  coagulation  case  (2  liters) 
the  concentration  is  2.5  *  10'  partic les/cc .  Je  the 
open  coagulation  case  the  total  volume  of  gases  pro¬ 
duced  by  the  7  g  of  propellant  is  calculated  to  be 
8  liters  at  300°K  for  a  seeding  particle  concentra¬ 
tion  of  6  x  1 06 pa r t ic  le s/cc . 

2.  Three  grams  of  BA  is  about  0.03  mole;  this  is 

1.8  x  1022  molecules,  giving  a  vapor  concentration 
in  the  confined  coagulation  case  of  1019  molecules/cc 
and  a  cone en t ra t i on  of  2.5  x  1018  molecules/cc  in  the 
open  coagulation  case. 


The  equation  for  the  collision  rate  between  molecules  of  two  dif¬ 
ferent  species  in  a  gas  is 


4  B 


-  V2  r mtuia  *  ,y~ 
4  L 


a ) 


whe  re 

Z, „  =  number  of  collisions,  sec  between  molecules  A  and  B 

A  o 

n  =  number  of  molecules /cr 

crA  B  ~  collision  cross  section  (diameter  of  molecules 
A  and  B) 

/?  -  gas  constant 
T  =  absolute  temperature 
M  ~  molecular  weight. 

If  Equation  (1)  is  applied  to  collisions  between  vapor  molecules 
and  aerosol  particles,  the  effect  of  a"v  assumed  molecular  weight  for 
the  average  aerosol  particle  is  essentially  canceled  by  the  form  of  the 
expression,  i.e.,  where  MA  «  Mg 


\  1 
M.  )  ~  "a 


(2) 
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The  calculated  values  of  the  collision  rates  between  BA  molecules  and 
silica  particles  are  ZAg  (confined)  *  2  *  1023  c ol 1 i s i ons/sec  at  500  K 
and  ZAS  (open)  =  1022  collisions/sec  at  300°K.  The  Smoluchowski  eouation 
expressed  as  a  collision  rate  is 

a  =  KnAnB  ' 

where  K  is  the  coagulation  constant.  In  systems  where  the  aerosol  par¬ 
ticles  are  all  the  same  size,  K  =  3  to  10  *  10_'°  cc/sec,  in  aerosol 
systems  containing  two  sizes  of  particles,  however,  the  K  for  collision 
between  the  two  sizes  increases  rapidly  with  the  ratio  of  the  sizes. 

Fuks1  has  calculated  values  of  K  for  various  size-pairs:  the  value  of 

K  for  O’  -  4  *  10-8  cm  (average  diameter  of  BA  molecule)  and 

a  =  2  x  10”4  cm  (average  diameter  of  silica  particle)  becomes 

KaS  -  2  *  10’4.  The  collision  rates  based  on  this  value  of  KAB  are  ZAB 

(confined)  =  5  x  1022  col 1  isions/sec ,  o..d  ZAg  (open)  =  3.0  *  1021  col- 
1  isions/sec . 

In  open  s l r g  1  e - s t age  coagulation  the  aerosol  cloud  is  rapidly  cooled 
and  the  critical  supersaturation  ratio  of  the  disseminated  vapor  is 
rapidly  exceeded.  This  produces  a  large  homogeneous  nucleation  rate, 
which  competes  wil.h  the  condensation  on  foreign  nuclei  and  produces  great 
numbers  of  very  small  (approximately  5  *  1 0- ‘  cm)  particles.  These  par¬ 
ticles  coagulate  rapidly  to  the  2  x  1 0” '  cm  size  and  become  visible  as 
a  persistent  smoke. 

In  the  two-stage  coagulator  the  minimum  temperature  is  controlled  by 
regulating  the  pressure  during  the  first  several  seconds  of  the  burning 
and  vaporization  process.  By  selecting  the  proper  release  pressure,  a 
minimum  temperature  can  be  reached  which  will  produce  a  saturation  ratio 
in  the  confined  vapor  just  exceeding  one,  thus  preventing  the  occurrence 
of  homogeneous  nucleation  of  the  vapor.  Also,  because  of  the  very  low 
supersaturation  value,  heterogeneous  nucleation  on  small  foreign  particles 
(t.e.,  carbon  nuclei)  will  be  inhibited  in  favor  of  heterogeneous  nuclea¬ 
tion  on  large  foreign  particles  (i.e.,  large  particles  are  present  in  the 
form  of  added  seeding  material). 
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Ill  H0M06ENE0US  NUCLEATION  STUDIES 


A.  Background 

The  pressure  of  a  bulk  liquid  in  equilibrium  with  its  vapor  is  pro¬ 
portional  to  the  absolute  temperature  T.  For  an  ideal  gas,  the  logarithm 
of  the  pressure  p  is  a  linear  function  of  l/T.  It  lias  been  observed  thaf 
in  a  vapor  phase  system  with  no  boundaries,  condensation  wi I !  not  neces¬ 
sarily  occur  when  the  temperature  decreases  to  equilibrium.  In  the  com¬ 
plete  absence  of  foreign  droplets,  ions  or  walls  the  supersaturation 
reaches  an  upper  limiting  value  beyond  which  condensation  will  always 
occur.  This  condensation  is  called  “homogeneous  condensation."  For  a 
terminal  temperature  Tj  the  critical  supe r sa t u ra 1 1  on  ratio  is  constant 
and  is  a  function  of  the  nature  of  the  compound.  The  degree  of  super- 
saturation  is  measured  by  the  saturation  ratio  5,  which  is  defined  as  the 
ratio  of  the  pressure  p  of  the  supersaturated  vapor  at  temperature  Tj  and 
the  pressure  pm  of  th<  saturated  vapor  at  the  same  temperature  T  j .  It 
has  also  been  observid  that  the  rate  of  condensation  is  small  until  5 
reaches  a  critical  value  after  which  condensation  heroines  so  rapid  that 
it  is  impossible  to  measure  the  condensation  rate  directly,  for  example, 
by  counting  the  number  of  drops  formed  per  unit  time  and  per  unit  volume. 
The  critical  value  of  5  at  Tj  is  rplled  the  critical  supr ' sa Lu r a L i on  ratio 
and  is  written: 


P  ® 


Earlier  theories  of  homogeneous  condensation  related  5  to  the  physi¬ 
cal  characteristics  of  the  compound.  These  theories  were  usually  based 
cn  the  experimental  observation  that  "the  controlling  process  in  conden¬ 
sation  is  the  nucieation  of  the  vapor  into  droplets."  To  explain  this 
observation,  Gibbs2  proposed  the  condition  which  must  exist  for  nuclei  to 
form  in  a  vapor.  Later,  the  rate  of  transition  from  a  vapor  to  liquid 
phase  was  studied  using  classic  rate  process  theory. 
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1.  The  Liquid  Drop  Concept 


This  concept  was  originally  introduced  by  Gibbs.2  He  proposed  that 
it  is  possible  to  supersaturate  a  vapor  because  the  formation  of  a  liquid 
drop  requires  not  only  the  work  needed  to  transform  the  vapor  into  a 
liquid  but  also  the  additional  work  required  to  form  a  stable  liquid  gas 
interface.  This  additional  work  acts  essentially  as  a  barrier  of  poten¬ 
tial  in  kinetic  processes,  and  it  slows  down  the  condensation  process. 

In  fact,  the  free  energy  of  formation  of  a  droplet  can  be  expressed  as 
the  sum  of  two  terms:  a  "bulk"  term  which  corresponds  to  the  free  energy 
of  formation  of  the  liquid  state  from  the  gaseous  phase  in  the  absence  of 
a  surface,  and  a  term  which  is,  the  free  energy  formation  of  the  surface. 

The  free  energy  difference  (per  mole)  is  related  to  5  by 

AG  =  RT  In  Pm  -  RT  In  p  -  -RT  In  S  (5) 


and  becomes  negative  for  p  >  p m. 

For  a  drop  of  radius  r  this  free  energy  difference  (per  drop)  is 
given  by- 


AG” 


(6) 


where  R  is  the  gas  constant  and  V0  is  the  molar  volume.  The  surface  term 
is  simply  written: 


A Gj  -  4TJr2a 


(7) 


where  cr  is  the  surface  tension. 

As  the  size  of  the  droplet  increases,  the  volume  term  increases  more 
rapidly  than  the  surface  term.  Con sequen tly,  a  critical  drop  size  will 
be  reached  after  which  any  increase  in  radius  will  cause  a  decrease  of 
free  energy,  and  therefore  spontaneous  condensation  will  occur.  It  is  in 
the  above  sense  that  the  surface  properties  act  as  a  bariier  of  potential. 
The  total  free  energy  of  formation  of  one  condensing  droplet  is  then 


£Tj° 


(8) 
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and  is  maximum  for  a  critical  value  r*  of  r  such  that 

(9) 


(10) 

This  value  AG  is  thus  the  barrier  which  must  be  exceeded  for  homogeneous 
nucleation  to  occur.  These  critical  size  droplets,  also  called  nuclei  or 
clusters,  are  present  in  supersaturated  vapors,  even  if  no  phase  transi¬ 
tion  occurs.  Note  that  Equation  (?)  is  the  well-known  Gibbs-Thomson 
relationship. 

As  already  mentioned  above,  Gibbs2  established  on  a  firm  basis  the 
equilibrium  conditions  that  must  exist  for  nuclei  to  form  in  a  vapor,  but 
he  did  not  study  the  conditions  for  transition  from  a  supersaturated  vapor 
with  a  few  nuclei  lo  a  condensed  vapor  with  a  very  large  number  of  drops. 
The  adaptation  of  rate  process  theories  to  the  problem  of  homogeneous 
nucleation  was  the  next  step  in  understanding  this  phenomenon.  This  is 
described  in  the  next  section. 


and 


AG 


1 6  ttq  3 
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2 •  The  Transition  Process 

The  first  attempt  to  describe  the  nucleation  rate  J  was  made  by  Volmer 
and  Weber.3  They  proposed  that  J  be  taken  as  a  product  of  two  terms — an 
energetic  and  a  kinetic  term.  The  energetic  term  n*  ,  defined  as  the  con¬ 
centration  of  nuclei  in  the  vapor,  can  be  related  to  AG*  by 

n*  =  r.le~^G'/kT  ,  (11) 


where  n,  is  the  number 
state  between  clusters 
frequency  of  collision 
duced  from  the  kinetic 


of  mol  ecu  1 es/cm3  (this  establishes  an  equilibrium 
of  different  size).  The  kinetic  term  ir  is  the 
cf  nuclei  per  surface  unit  477r*2  and  can  be  de¬ 
theory  of  gases: 


P 

w  =  - 

< 2  k  r  > 


(12) 
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where  p  is  the  supersaturated  vapor  pressure  and  m  is  the  mass  of  the 
nuclei.  Hence,  the  Volmer- Weber  equation  for  the  rate  of  nucleation  is: 


J 


(13' 


It  follows  therefore  that  n j  can  be  related  to  the  super. saturation  vapor p 
in  the  case  where  the  vapor  is  considered  as  a  perfect  gas: 


n 


1 


P 

kT 


(14) 


•  * 

and  that  A G  is  related  to  S  (compare  Equations  (5)  and  (9).  Because 
the  rate  of  condensation  becomes  very  rapid  when  S  reaches  5  ,  it  is  not 
possible  to  measure  J.  Consequently,  the  theory  assumes  instead  that  when 
one  drop  is  formed  per  cc  per  second  ( J  -  1),  S  has  reached  its  critical 
value,  and  for  this  condition  a  theoretical  value  of  S  is  calculated  and 
can  be  compared  with  the  experimental  value  of  S. 

Other  theories  have  been  proposed  which  differ  from  Volmer  and  Weber’s3 
initial  formulation.  These  theories  differ  mainly  in  the  expression  of  one 
or  both  of  the  two  terms — kinetic  and  energetic.  Volmer  and  Weber’s  theory 
was  refined  by  mu  1 1 i p 1 y i ug  tbs  frequency  tern  by  a  condensation  coefficient 
(d[  ^  1)  to  account  for  collision  efficiencies  (af  is  usually  taken  to 
equal  one ) . 

Becker  and  Doering*  refined  the  energetic  term  by  taking  the  steady 
state  instead  of  the  equilibrium  state  to  calculate  the  concentration  of 
nuclei  in  the  vapor.  They  considered  the  set  of  bimolecular  reactions 

Al  *  4,  ~ A2  (IS) 

A2  *  4,  ^==±  43  (16) 

A  ...  *  A.  =±  4*  (17) 

I  1  1  >  i 

where  4(  is  the  critical  nucleus  containing  i  molecules.  The  concentra¬ 
tion  n*  is  then  calculated  for  this  steady  state  and  is  written  in  the  form 

n *  =  Znle~*c',kT  (18) 
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Here  Z,  usually  called  the  nonequi li br i um  factor,  is  given  by 


The  size  of  the  critical  cluster  is  then: 

2  a 

5*  ) 3 


(24) 


More  recently,  Lothe  and  Pound6  proposed  that  the  free  energy  of 
formation  AG  should  be  corrected  by  a  term  AGj  '  to  account  for  the  motion 
of  the  droplet  as  a  whole  (AG  being  computed  for  a  critical  cluster  at 
rest).  This  term  can  be  calculated  from  the  rotational  and  transitional 
partition  functions  of  the  cluster.  Lothe  and  Pound  included  another 
correction  term  A G^'  arising  from  the  requirement  of  conservation  of  the 
degrees  of  freedom  of  the  cluster  from  the  vapor  phase  to  the  liquid 
phase.  This  term6  has  been  shown  to  be  approximately  equal  to  Ts,  where  s 
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is  the  molecular  entropy  of  the  liquid.  The  rate  of  nucleation  after 
correct  ieti  is  then  written  in  the  form: 


J  '  a  ZT(  4rrr*  2  ) 


-ac*  *r 


(  277**7-)  •'* 


where  P  =  e"AC  /tT  is  the  Lothe-Pound  factor.  An  error  of  4?f  was  found 
in  chocking  the  derivation  of  this  equation;  the  corrected  equation  is: 

r  .  ,2., 

P  V 1  4  nfi6 


which  reduces  to 


Here  A  is  the  Planck  constant  h  divided  by  2^7,  and  1  is  the  kinetic  moment 
of  the  critical  cluster,  equal  for  a  spherical  dr op  let  to 


5  ‘  -V  r 


Taking  J  =  1  and  Ol  r  *  1,  we  have  utilized  the  Lothe-Pound  equation 
[see  Equation  (25)]  a0  as  a  function  of  d,  a,  p^,  and  5  .  introducing  a 
new  function  \  independent  of  S  ,  such  that 


K  -  In  S  -  12  In  (In  S  ) 


where 


K  =  In 


24/2  k* 


k*  S 

In  -  *4  In  a  -  —  -  In  pwd  <  4  in  VT 


“  6.21132  ♦  4  In  a  -  —  -In  pwd  *  4  In  MT  (30) 

R 
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0 


( 31  > 


andX  »  l  n  S  yields 

X 3  *  X 2  <|t  ♦  \  -  —  In  a  -  1 2  1  n  X  f  —  In  (a  *  X  »  X  (l?  In  X  -  X )  J 

This  simplified  form  of  hquation  {2  b)  is  equivalent  to  the  l.o  t  h  e  -  Pou  n  d 
nucleation  rate  equation.  Solving  this  equation  for  specific  values  of 
a,  b,  and  X.  we  obtain  from  the  root  X  =  In  S  the  theoretical  value  of 
the  critical  supe i so cu ra t i on  ratio,  which  will  be  compared  to  the  measured 
values  of  S'  .  t*  may  be  calculated  by 


(32) 


To  measure  the  critical  supersa Lu rat l on  ratio,  the  following  techniques 
have  been  used:  the  expansion  chamber  technique,  the  diffusion  technique, 

and  the  jet  technique.  The  basis  of  these  techniques  is  to  cool  the  vapor 
very  rapidly  b>  enforcing  sudden  changes  of  temperature.  It  is  hoped  that 
in  this  way  the  vapor  will  become  supersaturated  before  it.  can  diffuse  and 
condense  on  the  wall  of  the  measuring  container  oi  chamber 


a.  Expansion  Chamber  Technique 

The  measurements  of  Wilson  and  Powell,  using  expansion  chamber 
techniques,  give  consistently  high  values  lor  the  critical  snpersatui ation 
ratios  of  water.  The  expansion  chamber  technique  consists  oi  cooling  the 
vapor  bv  suddenly  expanding  a  mixture  of  the  vapor  with  an  inert  carrier 
gas.  If  the  expansion  is  assumed  to  be  adiabatic  and  reversible,  the 
final  temperature  can  be  computed  from  the  ratio  of  the  volumes  before  and 
after  expansion,  from  the  initial  temperature,  and  from  the  heat  capacities 
at  constant  pressure  and  constant  volume,  respectively.  Unfortunately, 
under  conditions  of  reversibility  and  adiabarity,  the  expansion  ratio  can¬ 
not  be  higher  than  1.3  or  14.  Th i s  limits  the  range  of  temperature  that 

can  be  investigated  and  if  the  critical  supersaturation  ratio  is  too  large 

• 

(as  in  the  case  for  CC14,  for  example)  it  may  be  impossible  to  measure  -S 
by  this  technique.  Wilson  and  Powell8  measured  the  critical  supe rsa l u ra¬ 
tion  ratio  of  water  for  temperatures  ranging  from  260  to  3253K.  Critical 
supersaturation  ratios  of  some  organic  compounds,  such  as  n i t rorac t h an e , 
e t h y 1 ac e t a t e ,  and  a  few  alcohols,  have  b“en  measured  by  Volmcr  and  Flood.9 
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Clarke  and  Podehush  measured  5  for  water  and  benzene  using  helium  as  a 
carrier  /cas.  Their  reported  »elurs  of  5  a:e  small — probably  because  of 
li *=•  t  e  rogeneous  condetisat  ion  . 

h-  Diffusion  Technique  end  Measurements 

Recently  Katz  and  Ostermier"  have  used  this  technique  for  measur¬ 
ing  critical  supersaturat  i  on  ratios.  Their  results  are  in  agreement  v-  i  t  h 
those  of  Volmer  and  F 1 oou . 9  The  aiilusion  technique  was  originally  used 
to  determine  only  the  6’  for  solidification  of  a  supercooled  liquid.  This 
technique,  which  uses  the  actual  steady-state  condition  during  the  measure- 
merit,  appears  to  be  very  promising  for  obtaining  critical  supersaturation 
ratios.  Unfortunately,  it  is  also  limited  to  a  small  tempeiature  range, 
and  it  requires  the  knowledge  of  the  transport  coefficients. 

c.  The  Jet  Technique 

This  technique  has  been  used  by  Higuchi  and  O'Konski.12  It  is  a 
valuable  approach  to  critical  super  sat oration  measurement  but  it  can  be 
used  cnly  for  compounds  of  low  volatility  because  of  experimental 
l ini  tat ions. 

3.  Discussion  of  Comparative  Data 

The  critical  supersaturation  ratios  measured  by  Volmer  and  Flood9 
for  organic  compounds  are  found  to  be  in  excellent  agreement  with  the 
values  predicted  by  the  Becker- Doeri ng  equation.  The  Wi 1  son - Powe 1 1  mea¬ 
surements  for  water  were  also  in  a  relatively  good  agreeme  n  t  d  l  II  l  gf  fl 
t  perature.s. 

It  was  not  until  recently,  when  Lothe  and  Pound6  suggested  their  cor¬ 
rect,  ion  to  the  Becker- Doering  theory  (a  correction  that  increases  the 
niicleation  predicted  by  the  Becke  r- Doe  r  i  ng  equation  by  a  factor  of  1  0  1  7  )  , 
that  a  possibl*  means  of  reconciling  Becker- Doeri ng  theory  with  experiment 
was  seriously  attempted. 

In  view  of  these  discrepancies  it  became  important  to  plan  further  ex¬ 
perimental  investigations  to  determine  wl. ether  the  correction  proposed  by 
Lothe  and  Pound  was  necessary.  For  example,  it  appeared  desirable  to  mea¬ 
sure  .V  for  'inter  over  a  larger  temperature  range.  In  the  course  of  our  con¬ 
densation  studies,  we de v c 1  oped  a  me t ho J  of  measuring  5  over  very  large  temp¬ 
erature  ranges.  Thi  s de vciopccnt  and  the  5*  dat a acqui red  will  be  described 
in  the  following  sections. 
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B.  Experimental  Studies 


In  devising  a  new  technique  for  measuring  S‘  a  flow  method  was 
chosen,  which  was  similar  to  the  1-flow  method  used  in  gas  kinetics  to 
study  reaction  rates.  The  T-flow  method  has  been  described  by  Bodenstein 
and  Wnlgast,L'  Langmuir, w  und  Benton  :C| 

1.  Description  of  the  New  Apparatus 

Instead  of  mixing  the  flows  arriving  from  the  br«iuiies  oi  the  Y,  it 
was  decided  to  have  the  gases  arrive  head-on  at  the  mixing  zone.  The 
basic  principle  of  the  apparatus  is  simple.  A  hot  gas  stream,  saturated 
with  the  vapor  of  interest,  is  mode  to  collide  head-on  with  another  gas 
stream  which  is  much  colder.  The  two  streams  come  together  in  a  small 
chamber;  the  mixing  region  between  the  tube  ends  is  centered  in  the  path 
of  a  strong,  collimated  light  beam  which  is  placed  at  such  an  angle  as 
to  produce  forward  scattering  of  maximum  intensity  from  any  droplets 
which  might  be  formed.  The  cold  stream  temperature  is  controlled  by¬ 
passing  the  diluent  stream  through  a  massive  cylindrical  copper  thermo¬ 
stat  submerged  in  a  cooling  bath.  The  hot  copper  thermostat,  which  is 
heated  with  electric  heating  tape,  is  placed  directly  over  the  cold  one 
and  the  two  are  mutually  insulated  by  3/4-in.  cork.  The  hot  copper 
thermostat  is  also  an  isothermal  vaporizer.  The  liquid  of  interest  is 
placed  in  a  well  in  the  hot  cylinder.  A  measured  flow  of  nuclei-free, 
dry  air  is  passed  into  the  thermostat ed  liquid  through  a  fritted  glass 
bubbler.  The  gas,  vapor,  and  entrained  droplets  then  pass  through  two 
layers  of  Gelman  Type  A  glass  fiber  filter  to  remove  all  of  the  liquid 
droplets  and  tc  assure  saturation  of  the  gas  stream.  The  gas- saturated 
vapor  stream  is  then  passed  through  a  deionizer  in  the  form  of  a  small 
annular  condenser.  The  copper  inlet  line  to  the  cold  thermostat  has  a 
similar  condenser.  The  distance  between  the  surfaces  in  these  condensers 
is  1  n~,  and  they  are  about  5  cm  iong.  A  45- volt  field  is  used  for  gas- 
ion  removal.  Figure  3  shows  a  detailed  schematic  of  the  apparatus. 

2 .  Technique  Description 

Both  the  hot  and  the  cold  flows  must  be  limited  to  laminar  flows. 
This  provides  for  steady- state  conditions  in  the  mixing  region,  isolated 
from  any  wall  effects,  as  shown  in  Fig.  4.  Since  stream  temperatures, 
pressures,  und  composition  are  known,  only  the  heat  capacities  of  the 
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where  is  the  average  heat  capacity  of  the  air,  is  the  average 

heat  capacity  of  the  compound  vapor,  F  and  F a  are  the  cold  and  warm 
stream  flow  rates  and  T  and  Ta  are  the  cold  and  warm  stream  tempera¬ 
tures.  By  combining  Equations  (34)  and  (35),  Tf  can  be  calculated. 

When  7’^  is  established  S*  is  easily  calculated: 

(measured  p  at  T. ) 

S*  =  - ■ —  (36) 

(saturated  pm  at  Tf ) 

where  the  measured  vapor  pressure  in  mm  Hg  at  Tf  is 


(37) 


where  p'  is  the  saturated  vapor  pressure  at  Ta  and  pa  at  Tf  is  deduced 
from  the  curve  for  saturated  vapor  pressure  versus  1/T  for  the  pure  compound. 

To  gain  a  better  understanding  of  the  capabilities  and  limitations 
of  the  technique,  it  appeared  desirable  to  study  the  variation  of  5* 
with  the  ratio  of  temperature  Ta/Tf  or  the  ratio  of  flow  rate  Fa/Fc  of 
the  gases  being  mixed.  From  Equations  (33)  -  (37),  we  can  deduce  S  and 
T^,  and  we  notice  that  if  T f  and  Ta  are  kept  constant,  *S  will  vary  with 
the  ratio  F e/F  .  There  is  a  maximum  value  for  *S  that  will  correspond 
to  a  specific  value  of  Fc/Fa  and,  since  T(  and  Ta  are  kept  constant,  to 
a  specific  final  temperature  Tf .  Figure  5  shows  a  curve  of  -S  versus  1/T 
for  CCLj  at  Tr  =  175°K  and  T.  =  300°K.  Curves  (such  as  the  one  in 
Fig.  5)  can  be  used  to  visualize  how  5  could  vary  to  reach  the  observed 
critical  values,  and  to  choose  the  appropriate  values  for  Tf  and  Ta  for 
a  specific  experiment.  Thus,  it  can  be  seen  that  5  will  reach  this 
critical  value  by  varying  Fc/Fa  and  keeping  T(  and  Tp  constant  (see 
Fig.  5). 

The  interrelationships  between  accuracy,  -S* ,  and  the  F  c/F  a  ratio 
were  quite  complicated  because  of  the  physical  characteristics  cf  the 
system.  The  5*  is  easy  to  measure  when  an  initially  large  f  /F  ratio 
is  decreased  slowly  to  the  onset  of  condensation.  The  drawback  to  this 
approach  is  that  5*  occurs  at  large  fe/Fa  ratios  (cf.,  Ff /F  =  70  in 
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FIG.  5  EXAMPLE  OF  LIMlTb  OF  COLD  AND  WARM  FLOW  RATIOS 
IN  THE  MIXING  STREAMS  TECHNIQUE 

big.  5),  under  these  cond  i  t  i  oiis  ,  the  mixing  region  is  displaced  toward 

the  in.  aii  millrl.  and  this  may  al  led  l  lie  observed  condition  lor  the 

onset  of  condensation.  When  .S*  is  approached  from  the  other  direction, 
i  .  e  .  ,  when  an  initially  small  p '  /  F  *  ratio  is  increased  slowly  l.o  t  lie 
onset  ol  condensation,  the  value  of  the  F  /F  ratio  is  close  to  one 

C  * 

( cf F,Fm  1.0  in  Kig.  5).  and  the  mixing  region  is  not  displaced. 
However,  large  amounts  of  condensate  form  on  the  viewing  chamber  sur¬ 
faces  and  cause  d i  I  I  t  ••  u I t i es .  We  have  been  making  Measurements  using 

the  initially  small,  increasing  F  /F  approach  for  comparative  purposes 
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In  measuring  critical  supersaturation  ratios  when  Fc  =  Fa,  a  large 
temperature  difference  between  the  two  air  flows  is  necessary.  We  found 
that  for  Fc  *  F  ,  the  final  temperature  range  ia  limited  to  a  practical 
interval  of  about  60°K  (192  to  253°K  in  the  case  of  CCl^.  Thus  it 
remains  necessary  to  accept  the  experimental  difficulties  accompanying 
the  iarge  F(/Ff  ratio  approach  in  order  to  conveniently  measure  S*  at 
temperatures  outside  this  temperature  range. 

Because  of  the  relatively  small  dimensions  of  the  mixing  zone,  it 
is  net essary  to  determine  whether  or  rot  the  residence  time  of  the 
mixing  gas  and  vapor  was  adequate  to  establish  steady-state  conditions 
of  cluster  and  embryo  populations.  The  induction  period  (time  required 
to  form  clusters  and  embryos)  at  critical  supersaturation  is  estimated 
to  be  about  10“ 1  sec.  At  a  total  air  flow  of  2  liters/min.,  the  gas  in 
the  mixing  zone  travels  about  5  x  10"5  cm  in  10“7  sec.  Since  the  minimum 
dimension  of  the  mixing  is  its  thickness,  estimated  to  be  about  10“3  cm, 
the  residence  time  is  more  then  adequate  to  establish  steady-state 
cond i t i ons . 

Finally,  an  analysis  was  made  of  the  effect  of  both  total  air  flow 
and  the  ratio  of  the  cold  air  flow  to  the  warm  air  flow  on  the  accuracy 
of  critical  supersaturation  ratio  measurements. 

The  optimum  total  air  flow  should  be  between  1  and  2  liters/min. 

Flo  3  W  h  ich  are  too  small  cause  errois  in  measurement  of  the  temperature. 
Flows  greater  than  2  liters/min  produce  turbulence  in  both  the  cold  and 
warm  air  streams  and  could  cause  distortion  of  the  mixing  zone. 

Using  this  technique  S*  values  were  measured  for  benzene,  CCl^,  CS,, 
CHC1 3 ,  and  HjO.  It  was  possible  to  obtain  very  high  values  for  S*  for  a 
range  of  temperature  much  larger  than  any  range  that  had  been  investi¬ 
gated  in  previous  studies. 

3-  Potential  Sources  of  Errors  and  Description  of  the  Second  and 
Third  Generation  Apparatus 

As  measurements  progressed,  it  was  possible  to  gain  a  better  under¬ 
standing  of  potential  sources  of  errors  in  the  new  technique.  These 
potential  errors,  together  with  the  corrective  steps  that  were  adopted 
to  improve  the  apparatus,  are  discussed  below. 
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A  possible  source  of  error  in  the  very  cool  temperature  experiments 
could  have  been  created  by  water  from  the  ambient  wet  air  diffusing  into 
the  mixing  chamber  and  condensing.  To  eliminate  this  possibility  the 
apparatus  was  modified  to  isolate  the  condensation  chamber.  Thi'.  change 
eliminated  the  possibility  of  partial  obscuration  of  the  condensation 
tone  by  external  moisture  condensation.  The  modified  apparatus  is  shown 
in  Fig.  6.  The  first  condensing  droplets  were  much  easier  to  detect 
than  they  had  been  previously,  and  the  accuracy  of  the  apparatus  was 
improved  considerably. 

The  arrival  of  the  dry  air  directly,  without  its  being  brought  to 
the  equilibrium  temperature  of  the  upper  thermostat,  created  problems 
in  measuring  5*  of  compounds  with  low  or  high  vapor  pressure,  since  the 
temperature  of  the  upper  thermostat  was  quite  different  from  room 
temperature.  To  eliminate  this  limitation,  a  new  vapor  saturated- warm 
air  thermostat  apparatus  was  designed  and  fabricated.  The  general 
features  were  the  same  as  in  the  first  apparatus,  but  the  gas  preheating, 
the  thermos  tatir.g,  and  the  filtering  were  improved.  Figure  7  shows  the 
new  parts  in  detail.  The  air  flow  temperature  can  reach  the  thermostat 
temperature  or  the  equilibrium  temperature  of  the  saturated  vapor  before 
it  mixes  with  the  liquid  to  carry  the  vapor.  In  Fig.  7  the  air  is 
preheated  at  the  cylinder  temperature  by  circulating  it  in  a  copper  pipe 
around  the  copper  cylinder.  Then  it  goes  into  the  vaporization  chamber 
to  carry  the  satuiated  vapor  through  the  glass  tube  and  the  fritted 
glass  bubbler.  In  the  first  prototype  the  air  was  channeled  directly 
to  the  glass  tube. 

We  have  also  made  some  measurements  using  two  warm  sections  of  the 
apparatus,  set  up  in  series.  The  first  warm  section  was  used  as  a 
vaporizer,  producing  saturated  vapor  at  temperature  T  •  The  second  warm 
section,  into  which  the  saturated  vapor  was  introduced,  was  used  as  » 
thermostated  superheater.  This  arrangement  provides  a  simple  way  to 
clear  the  viewing  chamber  of  condensate.  This  can  be  done  by  circulating 
dry,  clean  air  in  the  system  before  admitting  the  vapor  to  the  epparatus. 
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C.  Experimental  hesults 


Curves  of  In  S*  versus  1/T  were  constructed  for  several  simple  com¬ 
pounds  over  wide  temperature  ranges. 

1  .  Wat  e  r 

The  critical  supersaturation  data  for  water  vapor  from  e  previous 
study16  and  from  the  present  work  are  summarized  in  Fig.  8-  A  comparison 
of  our  results  with  the  results  of  others  shows  that  there  is  agreement, 
between  our  measurements  and  those  of  Wilson  and  Powell,8  whose  results 
have  been  considered  the  best  available. 

Sander  and  Damkoehler17  and  Madonra  et  a! . ,  18  claimed  to  have  observed 
freezing  nucleation  of  water  at  lo’v  ten.f  features  (homogeneous  condensa¬ 
tion  from  the  vapor  phase  to  the  solid  state).  We  have  not  observed 
scintillations  upon  condensation  and  construe  this  to  mean  that  the  con¬ 
densed  particles  were  liquid.  As  the  temperature  decreased,  the  size  of 
the  nuclei  kept  getting  smaller,  and  the  aspect  of  the  cloud  kept  chang¬ 
ing  gradually  from  droplets  to  smoke.  The  temperature,  -65°C  is  generally 
believed  to  be  a  transition  tempeiature,  beiow  which  water  must  crystallize 


FIG.  8  COMPARATIVE  CRITICAL  SUPERSATURATION  CURVES  FOR  WATER  VA^OR 
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when  con  den  s  l  ng .  Hirth  and  Pound14  Lake  -41°C  as  this  lower  L  einpe  ra  Lur  e 
limit  because  -31°C  is  the  critical  temperature  for  rapid  homogeneous 
nucleation  of  ice  from  supercooled  water,  but  our  observations  lead  us 
to  believe  that  nucleation  occurs  from  the  vapor  to  the  liquid  state  at 
temperatures  lower  than  -65°C. 

2.  Organic  Compounds 

Figures  9-12  show  curves  of  In  5*  versus  1/7'  for  benzene,  carbon 
tetrachloride,  carbon  disulfide,  and  chloroform.  Computations  from  the 
experimental  data  were  done  with  the  TWX  Time  Share  Computer  system. 

Re  suits  of  an  error  analysis  shewed  that  both  saturated  vapor  pres¬ 
sure  and  temperature  had  been  measured  with  a  high  degree  of  accuracy  and 
that  largest  errors  were  probably  caused  by  promoters  of  nucleation  such 
as  foreign  nuclei.  Nevertheless,  we  did  find  some  systematic  error  in  a 
small  range  of  temperature,  and  we  investigated  the  possibility  that  it 
ight  have  resulted  from  a  temperature  gradient  in  the  cold  thermostat 
when  liquid  nitrogen  was  used  as  cooling  bath  Because  the  error  was 
systematic  we  were  able  to  correct  the  measured  T  (when  the  liquid 
nitrogen  level  was  beneath  the  copper  thermostat  a  maxinum  correction 
of  -20CK  had  to  be  applied). 

For  benzene  we  were  able  to  make  a  very  difficult  measurement  at  rela¬ 
tively  high  temperature  T  -  270°K  for  equal  to  /■'  approximately  and  found 
5*  -  8.0.  We  obtained  a  slower  rate  of  nucleation  than  Clarke  and  Hodel  sh,to 
who  found  S’  =  8  at  250JK.  This  result  indicates  that  they  may  have  h 
heterogeneous  condensation  when  they  measured  .S’  for  benzene.  The\ 
using  helium  as  a  carrier  gas. 

D.  Discussion  and  Conclusions 

1 .  Comparison  of  Fxperimental  Results  with 
Earlier  Theories 

Figures  9-13  show  the  curves  of  In  S’  as  a  function  of  the  inverse 
of  the  absolute  temperature  1 ff  for  water,  benzene,  carbon  Let rac h lot ide , 
carbon  disulfide,  and  chloroform.  To  calculate  the  theoretical  critical 
supersaturation  ratio  for  t.l  jse  compounds,  we  have  extrapolated  the 
1 i qu i d - 3 t a t e  physical  properties  to  very  low  temperatures.  The  experi¬ 
mental  data  appealing  in  the  theoretical  expressions  s't  the  density. 
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FIC.  10  LOGARITHM  OF  SUPER  SATURATION  RATIO  AS  A  FUNCTION 
OF  THF  INVERSE  OF  THE  ABSOLUTE  TEMPERATURE 
FOR  HOMOGENEOUS  NUd  EAT  ION  OF  CARBON  TETRACHLORIDE 


FIG.  11  LOGARITHM  OF  SUPERSATURATION  RATIO  AS  A  FUNCTION 
OF  THE  INVERSE  OF  THE  ABSOLUTE  TEMPERATURE 
FOR  HOMOGENEOUS  NUCLEATION  OF  CARBON  DISULFIDE 
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FIG.  12  LOGARITHM  OF  SUPERSATURATION  RATIO  AS  A  FUNCTION 
OF  THE  INVERSE  OF  THE  ABSOLUTE  TEMPERATURE 
FOR  HOMOGENEOUS  NUCLEATiON  OF  CHLOROFORM 
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the  saturated  vapor  pressure,  the  surface  tension,  and  the  molar  entropy 
at  temperature  T- 

The  logarithm  of  the  saturated  vapor  pressure  p ^  is  extrapolated 
graphically  from  the  linear  relationship  with  1/T  (see  Appendix).  The 
density  is  determined  by  the  equation 

d  «  otj  +  a 2t  +  a3t2  *  a4t3  ,  (38) 

where  t  is  the  temperature  in  degrees  centigrade  and  a  j  .  a2,  (t  ,  and  a  4 
are  experimental  constants  (see  Appendix).  For  c  at  T  we  use  the  equation: 

a  =  ^1  *  @2t  *  ^3f2  +  ’  139) 

where  ,  yfl  and  /3^  are  constants  determined  experimentally  (see 

Appendix).  The  liquid  entropy  s  is  computed  from 

C 

s  =  s°  *  J  — -  dT  ,  (40) 

T 

298.2  ' 

where  s°  is  the  standard  molar  entropy  for  liquid  state  (at  25°C  and 
atmospheric  pressure)  and  C  is  the  molar  heat  capacity  of  the  liquid  at 
constant  pressure. 

The  physical  property  data  were  compiled  for  the  organic  compounds 
and  water  and  then  used  m  the  nuclear  ion  rate  equations  of  Becker- 
Ooering  and  Lothe  and  Pound  to  obtain  the  calculated  values  of  In  S’, 

R  and  X  (see  Figs.  9-13). 

2(  In  S’)3  *  b( In  S’)2  -  a  «  0  ,  (21) 

X3  ♦  X2  <b  ♦  .V  -  -  In  a  -  12  In  X  ♦  -  In  [a  ♦  X3  ♦  X2  (12  In  X  -  A)]V  a  =  0 

XI  2  2  J  (31) 

A  comparison  of  the  theoretical  results  obtained  for  R  and  X  with 
the  experimental  values  of  In  S’  shows  a  discrepancy  between  theory  and 
experiment,  this  discrepancy  increases  very  rapidly  as  the  temperature 
decreases.  Furthermore,  the  valuta  of  S’  become  very  large  at  low  tem¬ 
peratures  and  calculated  cluster  si <e  becomes  very  small.  This  discrep¬ 
ancy  may  be  due  to  the  fact  that  macroscopic  properties  such  as  p,  d, 
and  a  are  assumed  to  be  applicable  to  very  small  clusters.  The  dep  nden  e 
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of  the  physical  properties  on  cluster  size  appears  to  be  very  important 
in  the  case  of  surface  tension. 


2.  Variation  of  Surface  Tension  with  Droplet  Curvature 

The  problem  of  the  variation  of  the  surface  tension  wit!;  dronJ.et 
curvature  has  been  studied  very  intensively  by  many  investigators,  hut 
the  relationship  between  surface  tension  and  curvature  has  never  been 
resolved.  Because  of  its  critical  importance  in  the  nucl'eation  model, 
a  study  was  made  of  the  surface  tension  of  clusters. 

We  determined  indirectly  the  surface  tension  (c  ,  .;s  a  function  of 
curvature  in  the  following  way.  We  introduced  a  parameter  p,  defined  as 
the  ratio  of  the  surface  tension,  corrected  for  the  curvature  (cr  )  and 
the  bulk  value  (cr,).  The  values  of  p  were  determined  by  replacing  the 
theoretical  values  of  S*  in  the  nucleation  rate  equation  with  the  mea¬ 
sured  values.  If  the  Lothe-Pound  theory  is  correct  (except  for  the 
value  of  cr)  ,  then  the  function  p{  i*)  -  or/cr  would  allow  the  theory  to 
predict  correctly  the  value  of  the  critical  supersaturation  ratio.  We 
thought  that  use  of  the  experimental  values  of  5*  and  of  the  theory  to 
compute  p(i*')  might  also  shed  some  light  on  the  dependence  of  surface 
tension  on  cluster  size. 

With  p  thus  defined  and  with  a,  b,  and  k  defined  by  Equations  (22), 
(23),  and  (30),  we  obtain 


(41) 


\r  -  A.  = 


a 

4  in  - 

a 


12  In  p 


where  ar ,  br ,  and  kr  are  the  values  of  a,  b,  and  k  for  cr  =  cr  or 

r 

°r  =  aP3  -  (44) 

,  r  1 

b  -  b  +  ~  In  p  (45) 
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K  *  12  In  p 


(46) 


Substituting  ar  ,  br ,  and  Kr  from  Equations  (44)-(46)  into  (21)  (for  the 
Becker-Doering  theory]  and  into  (31)  [for  the  Bee ke r- Doe r i ng  theory  cor¬ 
rected  for  the  Lothe  and  Pound  factor]  and  setting  E  =  2n  S* ,  we  obtain 


two  new  equations: 


E2  In  p  -  £2 ( b  *  2l)  =  0 


ap3  -  llE2  In  p  -  —  In  [ap3  -  ]  2  E 2  In  p  ♦  E2 (E  ♦  12  In  E  -  A.)] 

-E2  ^6  +  E  +  A-  —  /na_  12  In  £  |  “  0  (48) 

Equation  (47)  yields  the  values  of  p  which  we  shall  use  for  the  correction 
of  the  surface  ten  si  on  for  sm all  droplets  in  the  .  a  s  e  of  t  h ;  Be  ike i  - 
Doering  theory;  Equation  (48)  does  the  same  for  tht  Lothe  and  Pound  theor\. 

Since  we  wish  to  determine  if  the  discrepant)  beiween  thcoreticcl 
and  experimental  values  can  be  ascribed  to  the  use  ol  o t  ,  »e  also  com¬ 
puted  the  number  of  molecules  i*  in  the  critical  cluster  i  xamine  the 
function  *  i  i  * . 


The  results  of  our  compulation  for  p  and  t*  aie  give'  in  bigs.  14-18- 


3.  Cluster  Si  and  the  Su  i  f  n  c  c  Te  i  >  i  <»n  Cor  r  ec  t  i  on 

In  the  case  of  water  the  f nn<  t  l on  .  decrease*  with  the  site  of  the 
d  r  o  o  1  e  t  s  and  c.  on  s  en  uen  t  1  v  increases  w  t  )■  njrvsiMit  !«■!  and  Su  n  d  u  u  i  s  t 


have  developed  a  o a ng 1 i ng - bo nd  model  for  wulei  on  _  assumption  that  the 
effect  of  curvature  on  surface-fret  energy  is  due  »  he  bioken  hydiog 
bonding  of  the  surface  molecules.  The)  catcula  erection  of  sppi 

mutely  1.2  for  i*  =  1G  to  ilG-  Oui  iiumr  m  a  1 1 1  ted  resu  is  are  i 

qualitative  agieement  with  their  calculation,  e:  all)  foi  the  pcinl 

(S*  -  5.07,  T  -  276. 3K),  where  p  is  1.10.  l|owe»  foi  the  point 

(5*  6.3  *  1 0 4 ,  T  '  165.24K),  their  caltula ted  <  . r i ec  t  l  in  of  1.23  1 o i 

i*  =  10  does  not  agree  with  our  result  of  1 .73  lot  i*  =  9. 


•he  bioken  hydiogen 
rrection  of  sppiuxi- 
ted  reiu'ls  are  in 
all)  foi  the  ptml 
fo i  the  p o  i  u l 
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FIG.  16  PARAMETER  p  AS  A  FUNCTION  OF  i* 

FOR  THE  SYSTEM,  CARBON  DISULFIDE  +  AIR 


FIG.  17  PARAMETER  p  AS  A  r'JNCTION  OF  i* 
FOR  THE  SYSTEM,  CHLOROFORM  +  AIR 


FIG.  16  PARAMETER  p  AS  A  FUNCTION  OF  i* 
FOR  THE  SYSTEM,  WATER  +  AIR 


In  the  case  of  nonpolar  molecules  (benzene,  CCl4,  and  CS2  )  ,  we  find 
that  p  does  riot,  increase  with  curvature  and  that  p  is  approximately  con¬ 
stant  when  the  size  of  the  nuclei  is  large  enough;  but  for  a  smaller  size 
nuclei,  p  increases  markedly  with  curvature.  )n  the  case  of  CHClj,  for 
large  size  values  of  t*,  Chloroform  is  a  slightly  polar  compound  end  its 
P  -  i*  curve  resembles  the  non-pol sr  case  except  at  low  temperature  (small 
i*),  where  the  effect  of  polar  forces  may  be  evident. 

It  is  to  be  expected  thet  surface  free  energy  for  small  drops  should 
increase  with  any  or i cn t a t 1 oned  force  field  and  specifically  with  dipole 
moment.  In  the  case  of  benzene  we  may  ha_  an  o r i en t a t . ona 1  effect  lha 
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increases  p  with  curvature.  Benrene  is  a  fiat  molecule.  end  it  should 
be  more  difficult  to  get  the  surface  molecules  in  their  oriented  positions 
and  therefore  a  geometrical  hindrance  effect  will  strongly  increase  /'. 

This  will  explain  why,  although  the  benrene  has  no  (•'•clru-  dipole  moment, 
p  increases  rapidly  for  small  values  of  t*,  as  in  the  case  of  polar 
molecules. 
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GLOSSARY 


A 

A * 

l 

a 

c 

c  , 

p  V  i 

c 

p 

Cv 

d 

F 

c 

F 

V 

r 

ag 

AG° 

AG* 

ag*2' 

AG” 

AG*' 

AG*/ 

.  * 

i 

i 

j 

k 

K 

m 


M 


n 


A  nucleus  or  particle 

Critical  nucleus  containing  i  molecules 
Condensation  coefficient 

Average  heat  capacity  at  constant  pressure  (cal/moie)  of  vapor 
Average  heat  capacity  at  constant  pressure 
Average  heat  capacity  at  constant  volume 
Density  (g/cc) 

Cold  stream  flow  rate  (ml/min) 

Warm  stream  flow  rate 
Lothe-Pound  correction  factor 

Free  energy  of  formation  of  droplet  (per  mole) 

Free  energy  of  formation 

Free  energy  change  of  formation  of  critical  nucleus 
Free  energy  change  of  rotation  of  critical  nucleus 
Free  energy  of  formation  (per  drop) 

Free  energy  change  accompanying  rotational  and  translational 
degrees  of  freedom  in  critical  nucleus 

Free  energy  change  of  translations  of  critical  nucleus 

Planck  constant  divided  by  2rr 

Number  of  molecules  in  critical  nuclei 

Moment  of  inertia 

Nucleation  rate 

Boltzmann’s  constant 

Coagulation  coefficient 

Molecular  mass 

Molecular  weight 

Number  of  molecules/cc 
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n 

* 

n 

i 

ni 

N 

co 

P 

P  CD 

AQc 

A<2» 

r 

* 

r 

R 

P 

s 

S 

S* 

s° 

a 


Concentration  of  critical  nuclei 
Concentration  of  nuclei  containing  i  molecules 
Number  of  molecules/cm 3 
Avagadro’s  number 

Frequency  of  collision  of  nuclei  per  surface  unit 
Vapor  pressure  (mm  Hg) 

Saturated  vapor  pressure  (mm  Hg) 

Total  heat  gain  of  cold  stream 
Total  heat  loss  of  warm  stream 
Radius 

Radius  of  critical  nucleus 

Gas  constant  (Boltzmann  constant  x  Avagadro’s  number) 

Ratio  of  surface  tensions  cr6/crr 

Molecular  entropy 

Supersaturation  ratio 

Critical  super satur at  ion  ratio 

Standard  molecular  entropy 

Surface  tension  (dynes/cm) 

Surface  tension  of  drop  of  radius  r 
Eulk  surface  tension 

Collision  cross  section  (diameter  of  particles  A  and  B ) 

Absolute  temperature 

Initial  temperature 

Terminal  temperature 

Temperature  (°C) 

Molar  volume 
Nonequilibrium  factor 

Number  of  collisions/sec  between  A  and  B  particles 
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APPENDIX 


PHYSICAL  PROPERTIES  USED  IN  NUCLEATION  CALCULATIONS 

Density  (d)  Relative  to  Water  at  4°C 

Water:  20<b)  d  =  0.9994  +  1.896  x  10'4(t  +  4)  -  1.0417  x  10" 5 ( t  + 

-  6.95  x  10“ B( c  +  4 ) 3 

CC.1 4 :  21  d  =  1.63255  -  1.911  x  IO’3  t  -  0.690  x  10'6  t2 

Benzene:21  d  =  0.90005  -  1.0636  x  10“3  t  -  0.0376  x  10'6  t 2 

CS2 :  21  d  =  1.29272  -  1.481  x  IO'3  t  -  3.06  x  IO'7  t 2 

CHC1 3 :  20,21  d  =  1.52643  -  1.8563  x  io'3  t  -  0.5309  *  IO"6  t2 

-  8.81  x  10'9  t3 

where  t  is  the  temperature  in  degrees  centigrade 

Surface  Tension  (cr)  in  dynes/cm 

Water:20  cr  -  76.96  -  0  1 5  2  (  t  <•  8 )  +  0.7882  x  10~4(t  +  8)2 

-  2.56  x  10" 6 ( t  +  8 ) 3 

CC1  ^ :  20,21  cr  =  29.38  -  0.13975  T  *  3.875  *  10"  4 1  2 

Benzene: 20,21  <7  n  31.58  -  0.1  37  t  +  0,0001  t2 

CS2: 20,21  a-  =  35.28  -  0,15217  t  *  0.8333  *  10“4  t2 

CHC. 1 3 : 20,21  a  »  28.6  -  0 .  1 363 ( «  -  10) 

Liquid  Heat  Capacity  Cp  at  Constant  Pressure  in  cal/mole 

Water:22  Cp  =  218.765  -  2.0182  T  +  6.769  x  io-3  T2  -  7.575  x  10 
CC14:  23  Cp  =  219.18  -  2.2173  T  +  8.5583  *  1C3  T2  -  1.0833  >-  10’ 

Benzene:24  Cp  =  -  37.816  +  0,6104  T-  1 . 8938  *  10“ 3  T2  +  2.140  *  10' 

CS2:  25  cp  =  17.325  +  7.6  x  10~4  T 

CHC1 3 : 21,26  Cp  =  12.4822  +  0.05137  T 
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Standard  Entropy  (s°)  (at  25°C/atn  of  pressure  in  c  a_lL_perjnole  _and 


per  degree) 

Water: 

s  0 

-  16.73 

CCl*:23*26 

s° 

=  51.25 

Benzene : 24 

s° 

=  41.30 

CS2 :  25 

,  0 
o 

=  36.10 

CHC1 , :  26 

s° 

=  48.5 
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Homogeneous  nucleation  always  occurs  in  condensing  systems  of  high  vapor  con¬ 
centrations  even  in  the  presence  of  foreign  nuclei.  Critical  supnrsaturation 
ratios  of  a  number  of  compounds  were  measured  by  a  newly  developed  experi¬ 
mental  method.  It  was  demonstrated  that  the  classical  Beckor-Doering  theory 
is  inadequate  for  the  preparation  of  nucleation  models.  A  modified  model  is 
presented  showing  the  interrelationships  among  physics?,  properties,  chemical 
structure,  and  temperature . 
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